function f=velocity(epsilonR,muR)
% velocity(epsilonR,muR) - wave velocity of coaxial cable
%  epsilonR - relative dielectric constant of dialectric (1 for air)
%  muR - relative permeability of dielectric (1 for air)

c=2.998e8;
% c - velocity of light in free space

f= c ./ sqrt(muR .* epsilonR);
